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In the paper random-site percolation thresholds for simple cubic lattice with sites’ neigh¬ 
borhoods containing next-next-next-nearest neighbors (4NN) are evaluated with Monte Carlo 
simulations. A recently proposed algorithm with low sampling for percolation thresholds es¬ 
timation [Bastas et al, arXiv:1411.5834] is implemented for the studies of the top-bottom 
wrapping probability. The obtained percolation thresholds are pc(4NN) = 0.31160(12), 

pc(4NN+NN) = 0.15040(12), p c (4NN+2NN) = 0.15950(12), p c (4NN+3NN) = 0.20490(12), 
pc(4NN+2NN+NN) = 0.11440(12), p c (4NN+3NN+NN) = 0.11920(12), p c (4NN+3NN+2NN) = 
0.11330(12), pc(4NN+3NN+2NN+NN) = 0.10000(12), where 3NN, 2NN, NN stands for next-next- 
nearest neighbors, next-nearest neighbors, and nearest neighbors, respectively. As an SC lattice 
with 4NN neighbors may be mapped onto two independent interpenetrated SC lattices but with 
two times larger lattice constant the percolation threshold pc(4NN) is exactly equal to pc(NN). 

The simplified Bastas et al. method allows for reaching uncertainty of the percolation threshold 
value pc similar to those obtained with classical method but ten times faster. 
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INTRODUCTION 

Finding percolation thresholds pc and observing clus¬ 
ter properties near percolation threshold [1-4] are one 
of the most extensively studied problems in statistical 
physics. The beauty of percolation [5] lays both in its 
simplicity and possible practical applications. The lat¬ 
ter ranges from theoretical studies of geometrical model 
of the phase transition [6], via condensed matter physics 
[7], rheology [8], forest fires [9] to immunology [10] and 
quantum mechanics [11]. 

In random-site percolation model the nodes of lattice, 
graph or network are randomly occupied with a probabil¬ 
ity p. The critical probability pc separates two phases: 
for p > pc the system percolates, i.e. one may find a 
single cluster of occupied sites which extends to the bor¬ 
ders of the system; while for p < pc only smaller clusters 
exist. Usually, the finite size scaling theory [12-15] is 
employed for percolation threshold pc estimation. This 
requires checking properties of some quantity X(p,L) in 
the vicinity of phase transition as it depends on the linear 
system size L 

X{p-L) = L~ x -f((p-pc)L !/"), (1) 

where J-(-) is a scaling function, x is a scaling exponent 
and v is a critical exponent associated with the correla¬ 
tion length [1], Eq. (1) yields an efficient way for pc de¬ 
termination as L x ■ X{pc\ L) = .F(O) does not depend on 
the linear system size L. It means that curves L x -X{p ; L) 
plotted for various values of L should have one common 
point exactly at p = pc- Unfortunately, the results of 


computer simulations rather rarely reproduce a single 
common point of curves X{p\L ) unless the number N lun 
of prepared lattices is very high. 

Recently, Bastas et al. proposed efficient method for 
estimating scaling exponents x and percolation thresh¬ 
olds pc in percolation processes with low sampling [16, 
17]. According to Refs. [16, 17] instead of searching for 
the point where curves X (p; L) intercept each other one 
may wish to minimize the pairwise difference 

k{p\x) = [H{p-,Li) - H{p-,Lj)] 2 , (2) 

with respect to both parameters x and p , where 

H( P ',L) = Y(p;L) (3a) 

as suggested in Ref. [16] or 

H(p-L) = Y(p-L) + l/Y(p;L) (3b) 

as proposed in Ref. [17] and in both cases 

Y(p-L) = L x -X(p;L). (4) 

The minimum of A(p; x) is reached for p = pc and x = 
fl/v, where /3 is a critical exponent associated with the 
order parameter (for instance probability of an arbitrary 
site belonging to the infinite cluster [1]). 

In this paper we propose simplified version of Bastas 
et al. algorithm, where only a single-parameter function 
A (p) must be minimized in order to provide percolation 
threshold estimation. With such approach we estimate 
simple cubic (SC) random-site percolation thresholds for 
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FIG. 1. Single sites from various neighborhoods of an SC 
lattice. The full neighborhoods contain z = 6, 12, 8 and 6 
sites for NN, 2NN, 3NN and 4NN neighborhoods, respectively. 

eight complex neighborhoods containing next-next-next- 
nearest neighbors. Our results enhance those of the ear¬ 
lier studies regarding percolation thresholds for complex 
neighborhoods on square [18] or SC [19] lattices. 

APPROACH 

Our proposition is to apply Bastas et al. technique 
for such quantity X (p; L) which does not require scaling 
along X axis by a factor L x in order to achieve statis¬ 
tical invariance of the shape X (p; L ) for various values 
of L. An example of such quantity is the (top-bottom) 
wrapping probability [2C ]: 

W(p-,L)=N(p- : L)/N Iua , (5) 

where IV (p; L ) is a number of percolating lattices with 
p ■ L 3 occupied sites among N lun lattices constructed for 
fixed values p and L. In thermodynamic limit we have 
W(p < pc\L —y oo) = 0 and W(p > pciL -A oo) = 1 
and thus scaling exponent a; of IP is equal to zero [1], 
Consequently, instead of the form given by Eq. (1) wrap¬ 
ping probability obeys a simplified scaling relation [1,2 ] 

W(p-,L)=g^p-p c )L 1 / v y ( 6 ) 

Eq. (6) again makes possible pc determination as 

W(pc',L) = 1/(0) does not depend on the system size 

L. 

Now, the equivalent of Eq. (2) may be written as 

A (p) = Y / [H(p;L i )-H(p-L ] )}\ (7) 

where 

H(p-L) = W{p-,L) + l/W{p;L). ( 8 ) 

Following Bastas et al. technique one should minimize 
function A (p); the found minimum may be then used for 
the pc estimation. 


Several numerical techniques allows for clusters of 
connected sites identification [21-24]. Here we apply 
Hoshen-Kopelman algorithm [22], which allows for sites 
labeling in a such way, that occupied sites in the same 
cluster have assigned the same labels and different clus¬ 
ters have different labels associated with them. 

Here we investigate an SC lattice with sites’ neigh¬ 
bors ranging from the nearest neighbors (NN), via 
the next-nearest neighbors (2NN) and the next-next- 
next-nearest neighbors (3NN) to the next-next-next- 
next-nearest neighbors (4NN). A scheme showing 
only single sites of each of the neighborhood types 
mentioned above is presented in Fig. 1. The full 
neighborhoods contain z = 6, 12, 8 and 6 sites for 
NN, 2NN, 3NN and 4NN neighborhoods, respectively. 
Also all available combinations of these neighbor¬ 
hoods are considered, i.e. (4NN+NN), (4NN+2NN), 
(4NN+3NN), (4NN+2NN+NN), (4NN+3NN+NN), 

(4NN+3NN+2NN) and (4NN+3NN+2NN+NN) con¬ 
taining z = 12, 18, 14, 24, 20, 26 and 32 sites, 
respectively. 

RESULTS AND DISCUSSION 


TABLE I. The critical values of pc for various neighborhoods 
based on minimization of A (p) function. 


neighborhood 

z 

Pc 

4NN 

6 

0.31160(12) = pc(NN) 

4NN+NN 

12 

0.15040(12) 

4NN+2NN 

18 

0.15950(12) 

4NN+3NN 

14 

0.20490(12) 

4NN+2NN+NN 

24 

0.11440(12) 

4NN+3NN+NN 

20 

0.11920(12) 

4NN+3NN+2NN 

26 

0.11330(12) 

4NN+3NN+2NN+NN 

32 

0.10000(12) 


For each pair (p, L ) of parameters N lun = 10 4 lattices 
with randomly occupied p-L 3 sites were simulated for L = 
40, 80, 120 and 160. The wrapping probabilities W ( p ; L) 
for various neighborhoods combinations are presented in 
Fig. 2. 

As it was mentioned in the Introduction, the results 
of computer simulations rather rarely reproduce a sin¬ 
gle common point of curves W (p; L) unless the number 
7V r un of prepared lattices is very high. It means, that 
finding the common point of W (p; L ) curves for various 
linear system sizes L may be quite problematic. In or¬ 
der to illustrate better this situation, we plot W (p; L) 
dependencies near pc with sites occupation probability 
step Ap = 10 -4 (see Fig. 3). And indeed, except for 
4NN+2NN+NN neighborhood, the curves W (p; L) for 
various pairs of L intersect at different points. Moreover, 
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FIG. 2. (Color online). Wrapping probability W(p;L) vs. occupation probability p. The results are averaged over N run = 10 4 
runs. The symbols (+, x, □) indicate the system linear sizes (L = 40, 80, 120, 160), respectively. 


for the smallest values of L the dependencies W(p;L) 
do not even increase monotonically with p. At the same 
time, the dependencies A(p; L) prepared with the same 
accuracy A p = 10~ 4 exhibit single and sharp minimum 
(it is worth to mention, that values of A are presented 
with the use of logarithmic scale). The minimum of A (p) 
corresponds to the percolation threshold pc■ The esti¬ 
mated thresholds are presented in Tab. I. 

The plots W{p]L) presented in Fig. 3 make possi¬ 
ble to detrmine the length of interval where the true 
value of percolation threshold is located. This length 
is equal to Sw{pc ) = 0.0004. Assuming that a real 
percolation threshold value is uniformly distributed in 
this interval one may evaluate the percolation thresh¬ 
old uncertainty as uw(pc) = bw{.Pc)/V?> ~ 0.00023. 
The approach based on A (p; L) dependence gives evalua¬ 
tion of 5\{pc) twice smaller and consequently u\(pc) ~ 
0.00012. On the other hand the method of pc estimation 
based solely on W(p; L ) dependencies, applied for similar 
neighborhood geometries leads to twice smaller lengths 
S(pc) and consequently similar uncertainties u(pc) but 
for ten times larger sampling (N lun = 10 5 ) [19]. One 
can conclude, that the method used by Bastas et al. 
leads to reaching uncertainty of the percolation threshold 
value pc similar to those obtained with classical method 
u\{p c ) ~ uwipc) but ten times faster. 

Note, that an SC lattice with 4NN neighbors may 
be mapped onto two independent interpenetrated SC 
lattices but with two times larger lattice constant. 
Thus we expect the percolation threshold pc(4NN) for 
the next-next-next-next-nearest neighbors to be equal 
exactly to pc(NN). Indeed, the obtained value of 
pc(4NN) «0.31160 agrees very well with values of per¬ 
colation threshold estimated for the nearest neighbors 


pc(NN) ss 0.31160768(15) obtained very recently in ex¬ 
tensive Monte Carlo simulation [25] and its earlier esti¬ 
mations [26]. 

Note, however, that reaching such accuracy requires, 
for L < 128, sampling over N lun = 5 • 10 8 lattices re¬ 
alization [25], while we recovered the first five digits of 
pc(NN) with statistics lower by more than four orders of 
magnitude. 

Knowing percolation threshold may be practically use¬ 
ful for many systems with neighborhoods ranging be¬ 
yond nearest neighbors [27] or next-nearest neighbors 
[28]. Thus practical application of pc values for longer 
ranges of interaction among systems’ items cannot be 
generally excluded in all typical applications of the per¬ 
colation theory, i.e. physics, chemistry, biology and social 
sciences. 

The work was supported by the Polish Ministry of Sci¬ 
ence and Higher Education and its grants for scientific 
research and PL-Grid infrastructure. 
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FIG. 3. (Color online). Wrapping probability W{p\L) and the pairwise sum A (p) vs. occupation probability p. The results are 
averaged over N Tun = 10 4 runs. The symbols (+, x, *, □) indicate the system linear sizes (L = 40, 80, 120, 160), respectively. 
The minima of A (p) correspond to the percolation thresholds pc- 
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